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Abstract 



We formulate the description of the dynamic confinement of a single fer- 
' mion by a flux tube. The range of applicability of this model extends from 

the relativistic corrections of a slowly moving quark to the ultra-relativistic 
^ \ motion in a heavy-light meson. The reduced Salpeter equation, also known as 

5^ , the no-pair equation, provides the framework for our discussion. The Regge 

structure is that of a Nambu string with one end fixed. Numerical solutions 
are found giving very good fits to heavy-light meson masses. An Isgur-Wise 
function with a zero recoil slope of ■^'(1) — —1.26 is obtained. 



1 Introduction 



The relativistic flux tube (RFT) model's |l|] major success has been its ability to 
account for leading relativistic corrections in heavy quark bound states and also to 
give the Regge structure of the Nambu string. The "Wilson loop low velocity expan- 
sion" of Eichten and Feinberg provides a rigorous non-perturbative framework 
for the discussion of these corrections. The spin-independent relativistic corrections 
were computed in an extension of this program @. The spin-independent QCD 
relativistic corrections are accounted for in a simple djTiamical picture by the RFT 
model IQ. The spin-dependent long range corrections have also been considered in 
the RFT model [Q. In the present work we show that the spin-orbit interaction 
obtained as a relativistic correction, differs from the "Thomas term" conventionally 
associated with QCD 0. 

Yet to be completely explored is the ultra-relativistic sector of the RFT model 
with fermions, although the spinless quark results have been encouraging ^ |[. In 
fact, with reasonable assumptions the Wilson area law exactly reproduces the orig- 



inal spinless RFT model [§, |T0|. One of the attractive features of the RFT model 



is that ultra-relativistic quark motion has little influence on meson's rotational dy- 
namics. For large angular momenta the tube's momentum and energy dominate and 
the RFT reduces to the Nambu string. The RFT flux tube can also stretch radially 
which involves no change in tube's momentum. For s-waves it then resembles a po- 
tential model. If both radial and angular motions are present Coriolis force couples 
the radial and rotational motions and implies a straight leading Regge trajectory 
accompanied by parallel equally spaced daughter trajectories [jT], 0]. 

In recent work []rT| we have examined the motion of a fermion in a confining 
central potential field and found that the Dirac equation is appropriate for scalar 
confinement while the "no-pair" variant of the Dirac equation (or reduced Salpeter 
equation), must be used for a Lorentz vector confining potential. Since the RFT acts 
like vector confinement for low angular momenta the Salpeter equation is relevant. 
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Our formalism for fermionic quark confinement is unusual in that the confinement 
is introduced into the kinetic rather than in the usual interaction term. The flux 
tube contributes to both energy and momentum to the dynamics so it makes little 
sense to consider it as a "potential" type interaction. By a covariant substitution we 
add the tube to the quark momentum and energy. We may equivalently view this 
as a "minimal substitution" of a vector interaction field. The result nicely reduces 
to the Nambu string in the limit in which the quark moves ultra-relativistically. 

In Section ^ we discuss the quantized form of the tube operator and the flux 
tube transformation. Starting from the reduced Salpeter (or no-pair) equation, we 
develop the radial equations with tube confinement in Section ^. A semi-relativistic 
reduction is presented in Section ^. In Section ^ the Regge structure of the RFT 
model for an ultra-relativistic fermion is determined, yielding a straight trajectory 
with slope a' = In Section |^ we present our numerical results. We obtain an 
excellent fit to the observed particle spectrum and compute the Isgur-Wise function 
for the semileptonic B D^*Hui decays. Our conclusions are found in Section |^. 



2 Flux tube transformation 

From the Lagrangian formulation of the spinless RFT model |l| it follows that the 
canonical linear momentum is sum of quark and the tube parts 

P = Pg + Pt , (1) 

where 

P, = — = Wr^^V^ , (2) 

r 

Pt = ^ = 2arF{vj_) , (3) 

with 

^ , , 1 ,arcsint>_|_ 1 , , 
F{v^) = — - — , 4 

4:V± v± 7_L 
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and the tube energy is found to be 



arcsint;_L . . 

iff = ar . (5) 

v± 



In these expressions 7^ = , „ and Wr = \ pi + m^- 

In an earher work the flux tube with fermionic ends was introduced through 
the free Dirac Hamiltonian by the covariant transformation, 

P'^v'-V't, (6) 

where 



= {HuVt), (7) 

Pi = (-r X t)pt = -r X — . (8) 

r 

Here, unlike in 0, we want to quantize this classical transformation at the outset. 
The prescription we follow is to first quantize iff, pt, and Pq by symmetrization of 
the classical expressions , so that 

a , arcsin v± , 
= -,{r,^^}, (9) 

Pt = y = a{r,F{v^)}, (10) 
P, = ^ = l{Wr,7±v±} , (11) 

with {A, B} = AB + BA and = — ^^r. The quantum operators corresponding 
to Pt and Pq are then 

Pt = ^(— xf-fx— ) (12) 
2 r r 

Pg = rpr + \{— xr-rx ^) (13) 
2 r r 

where Pr = is hermitian. From p = Pq + Pt and the identity 



we obtain the usual 



r —t 



dr' 



„ L 
r X — . 

r 



(15) 



The full flux tube transformation (P) of the free Dirac Hamiltonian is then given by 



where 



Ho{p) = oc p + (3m 



(16) 



(17) 



We choose to work with the reduced Salpeter equation |T^, 13 1 in the limit where 
the heavy antiquark mass is infinite . In coordinate space it is given by 



A+[E,-//o(p)-V^(r)]A+^(r) = 



(18) 



and also known as the no-pair equation [|T5|. Here, Eg = M — mq is the energy of 
the light degrees of freedom, V{r) is the short range vector interaction, while A+ is 
the positive energy projection operator 



with 



A. 



En 



2Eo 



(19) 



ip^ + w?. (20) 

The fiux-tube transformation (|16D on equation ([T8|) results in the no-pair equation 
describing heavy-light mesons with the flux tube accounting for the long range in- 
teraction 

[Eg - Ho{p) - A+{V{r) + Ht-cx- Pi)A+]^(r) = . (21) 
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3 Radial equations 



Coupled radial eigenvalue equations follow from the spherical symmetry of equation 



IT). As in 111] the solutions are of the form 



^'^ (r) 



where 




tk 

9j 





iy^'' 

The Dirac quantum number k labels the bound state and is defined by |]16 



The identities 



lead to 



^ jm 



crvy 



jm 



-yr'' 

^ jm ' 



(22) 



(23) 



(24) 



(25) 
(26) 



where 



Tt 



CT ■ p3^S(r) 
<^ ■ Pt3^S(r) 



,k .d 1, 



k) 



T^ySi^)Tt , 



-k{l - k) 



Pt 



1 + ^) k 
■Pt 



k{l + k) 



(27) 
(28) 



(29) 
(30) 



The no-pair equation can now be expressed in radial form. Starting from (^) and 
the wave function (p^, we use the above relations to obtain matrix radial no-pair 
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equations 



with definitions 



J 3 

9j 



m 



-m 



L 



2E, 



2E2 



V{r) + 



Tt 

V{r) + H-^ 



In the above expressions 



Ah 



Ei±m 



^0 — v"^ ~ -^T-' 

From (^4|) we observe that the tube energy enters much hke the time component of 
the short range Lorentz vector potential V{r), which we take to be 



D^D± . 



(31) 

(32) 
(33) 
(34) 

(35) 
(36) 



V{r) = - 



K 



(37) 



Finally, the only unknown operators in (^) are f which are contained in pf^ and 
H^^. These will be determined from the heavy-light orbital angular momentum 
equation as in the spinless RFT model [§. From 

-L 1 



:{Wr,^±v^} + a{r,F{v^)} 



v^(r) 



(38) 



and using Ly^j^^r) = yk{k + l)yj^{r), one immediately finds equations that deter- 



mine the v^^ operators: 



^~fc(-fc+i) 



l{Wr,^ivi} + a{r,F{vi)} 



(39) 
(40) 
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The numerical technique used to solve for v_[_ is discussed in detail elsewhere [0, ^ . 



4 Semi-relativistic reduction 



For comparison with previous results we decouple (0) and make the non-relativistic 
reduction to order (f/c)^ (ignoring the short-range V(r) component), 



H = E„ 



m + 



2Et 



{E+ + m)H^{E+ + m) + {E^ + m)TtD, 
1 D^D^ 



+ 



D_TtiE+ + m) - D_H-^D+ 
1 



+ D_TtD_D+- D_H^\E, 
From and (|1^) we find (with L ^ mv±r) 

Ht 

Pt 

and from (^) it then follows 



2E^ E^ +m 
{E+ + m)H^D_D+ + {E+ + m)Tt{E^ - m)D^ 

1 



m)DA 



2E^{E^ + m) • 



1 

~ ar H — art>_L ~ ar H — 

6 Gm^'r 

1 aL 

~ -arv± ~ , 

3 3m 



m 



2m 



5m"' 



+ 



1 



4m^ 



2m 



Taking expectation values, and using 



k{k + l) - 
-{k + 1) - 
D^D+ -p 



2S-L , 

2 



(41) 



(42) 
(43) 



(44) 



(45) 
(46) 
(47) 
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we obtain the effective Hamiltonian 

H = m + - ^ + ar — . (48) 

2m Sm-^ Qm^ r Qm^ r 12m^r 

Physical interpretation can be estabhshed for each of the terms of the above effective 
Hamiltonian. The first three terms result from the expansion of the relativistic 
energy \/p^ + m^, and the fourth is the static tube energy. The spin-orbit terms is 
the net result of the kinematic Thomas precession and tube energy and momentum 
contributions. As already mentioned, our result differs from the one conventionally 
found from QCD §. 

The term in (|S|) has a nice interpretation as pointed out previously |^ . The 
tube's explicit contribution to the angular momentum is the sixth term of (]48|). 
There is also a hidden contribution from the /2m term. The rotational part of p 
yields to leading tube order {pt ^ Pq) 

2m 2m m r 

= ^ + (49) 
2m mr Sm 

where we have used L ~ p^r and Lt/r c:^ aL/3m from (H3|). The total tube contri- 
bution to the part of the Hamiltonian is then 

aL^ ,1 1, aL"^ 
m^'r 6 o bm^r 

The above is exactly the rotational energy of a rod of a mass ar and length r rotating 
with angular velocity u = L/mr"^, 

* 2 

= \\ary^ = ^ . (51) 
2 3 m^r'^ 6m^r 

The final term in (|48| ) is a Darwin like term and it is sensitive to how the Hamiltonian 

is symmetrized. If one symmetrizes at the end as in previous work 0, a different 

result is obtained. 
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5 Regge structure 



A confined particle in a sufficiently large angular momentum state moves ultra 
relativistically. For linear confinement we expect linear Regge trajectories from 
classical considerations @, [l^]- In the flux tube model the angular momentum and 
rotational energy of the quark are negligible compared to those of the tube. The 
slope of the Regge trajectory can be found from the model without a detailed solution 
of the wave equation. To determine the Regge slope directly from the no-pair tube 
equation we note from (H) and (^) that (m — > 0, — > 0, v± ^ 1) 

p1 (52) 

r 

Tt -kpt , (53) 

with k = -|^ being the sign of k, while from @ it follows 

^ . (54) 

We also have 

Ei p, , (55) 

A- i, (56) 



SO that in the limit \k\ ^ 1 we obtain 




^0 — ' ^* I A; ' ' ^^^^ 



(58) 
(59) 
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Multiplying out LIL then yields 

Independent of the sign of k, the eigenvalues of (|60D are 

= H,^^, (61) 

= -Pt^--^- (62) 

The positive energy solution has the proper Regge behavior. From (|38D we find the 
total orbital angular momentum of the system in the ultra-relativistic limit to be 

|fc|»i avrr^ 

L — ^ - \k\ . (63) 

Eliminating r and using (|HT|) we obtain the ratio of angular momentum to the square 
of the energy of the light degrees of freedom, 

This is the same slope that one would obtain from a Nambu string with one end 
fixed. 



6 Numerical results 

By expanding and g'^ in terms of a complete set of basis states, and then trun- 
cating the expansion to the first states we can transform the angular momentum 
equations (p9|-|40|) into two N x N matrix equations from which matrices for v^'^ 
operators are found Also, the radial no-pair equation becomes a 2N x 2N 
matrix equation, which we solve using the Galerkin variational method. The quasi- 
Coulombic basis states (which depend on a scale parameter (3), and all matrix ele- 
ments used are described in [|ll|. In Fig. we show the dependence on jS for the 
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three lowest positive and three highest negative energy states, with the particular 
choice m = 0.3 GeV, k = —1, j = ^, a = 0.2 GeV^, and k = 0.5. As we increase the 
number of states A^, the plateau where the positive eigenvalues are stable enlarges. 

The Regge slope of the no-pair equation with the flux tube was shown in Section 
|]to be a' = Our numerical solution agrees as shown in Fig. |^. In this figure we 
illustrate the leading trajectories for the two light degrees of freedom parity states 
corresponding to k = ±(j + |) for rriu^ = 0.3 GeV, a = 0.2 GeV"^ and k = 0.5. We 
also show several daughter trajectories corresponding to radial excitations. 

In order to recover the universal Regge slope, we fix a to be 0.2 GeV^, and we 
also choose rriu^d = 0.3 GeV and fit to the spin averaged heavy-light meson states. 
Our result is shown in Table and parameters of the fit are 

mu,d = 0.300 GeV (fixed) , 
rris = 0.581 GeV , 
rric = 1.286 GeV , 

rub = 4.621 GeV , (65) 
a = 0.200 GeV^ (fixed) , 
n = 0.641 . 

The agreement of the fitted levels to experiment is excellent, all spin-averaged states 
are reproduced with errors of about 5 MeV. 

Finally we use the wavefunctions with parameters of ( |65D to evaluate the Isgur- 
Wise function for the semileptonic B D^*Hi'i decays using [|1^ 



«.) = ^(,„(2£j^.)), (66) 



{A) = / dr r^R{r)A{r)R{r) . (67) 
Jo 



where 



As shown on the Fig. y, the agreement with ARGUS pO| and CLEO II data is 
excellent. For the seven CLEO II data points, this IW function has P^r degree 
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of freedom of 0.4 (corresponding to about 90% CL). To calculate the slope, we use 
ID| expression 

1 1 

2 "^3' 



= -{7^ + -,El{r')) ■ (68) 



For the light quark mass of 300 MeV , and using values for a and k from ( |65D we 
find 

= -1.26 . (69) 

We observe that slopes calculated from the different versions of the no-pair equa- 
tion (with vector confinement |Tl]| or with the flux tube) are more negative (^'(1) 
ranging from about -1.2 to -1.3) than the ones calculated from the spinless RFT 
model or the Dirac equation with scalar confinement (^'(1) — —0.9). 



7 Conclusions 

In this paper we further explore an idea ||^ for describing the relativistic quantized 
system of a flux tube with fermionic quarks at its ends. We work out in detail the 
simplest case where one fermion has infinite mass. Our technique for incorporating 
the tube into the reduced Salpeter equation is the covariant flux tube transformation 
([161) . In order to achieve a satisfactory description of a heavy-light system, there 
are two physical requirements that have to be satisfied: 

1. In the ultrarelativistic limit the tube dynamics must dominate giving a Nambu 
string limit. 

2. The semi-relativistic corrections must agree with rigorous QCD expectations. 

This model meets the above requirements with the possible exception of the spin- 
orbit interaction. 

Let us briefly review the evidence concerning the QCD spin-orbit interaction. 
The conventional wisdom is that the spin-orbit effective Hamiltonian is equivalent 



13 



to pure Thomas precession implying a coefficient — | instead of — | in the S ■ L term 
of (^Sp. The analysis of Eichten and Feinberg |^ does not say what the spin-orbit 
coefficient is although later work by Gromes p and the Milan group P], with further 
assumptions, find the Thomas result. We believe this remains an open question. 

The model also yields the expected Regge structure. For large angular momen- 
tum a Regge slope of a' = characteristic for a Nambu string with one fixed end, 
is obtained. A sequence of parallel equally spaced daughter trajectories follows from 
radial excitations of the quark. By varying the tension a, the Coulomb coefficient 
K and the quark masses we can account for the known spin-averaged heavy-light 
levels. In Table |1| we show a fit with the light quark mass and the tension fixed at 
reasonable values. We also used the s-wave wave function to compute the Isgur-Wise 
function. The result agrees well with the known experimental data, as shown in Fig. 

i 

The numerical techniques developed for the flux tube model with spinless quarks 
[]T|, ^ can be directly applied in the fermionic case since quantization of the orbital 
angular momentum sector is independent of the spin part. The extension to arbitrary 
two fermion systems seems to be free of any intrinsic difficulties. 
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TABLES 



Table 1: Heavy-light spin averaged states. Theoretical results are obtained from 
the no-pair equation with the flux tube. Spin-averaged masses are calculated in 
the usual way, by taking | (|) of the triplet and | (|) of the singlet mass for the 
s(p)-waves). 



state 


spectroscopic label spin-averaged 


q- 


n. 


theory 


error 






mass (MeV) 


j 


A; 


(MeV) 


(MeV) 


cu, cd quarks 
















D (1867) 
D* (2009) 


0- 

1- 




IS (1974) 


1 

2 


— i 


iyoU 


D 


Di (2425) 
D* (2459) 


1+ 
2+ 


'Pi' 

'P2 


IP (2446) 


3 
2 


-2 


2440 


-6 


cs quarks 
















Ds (1969) 
D* (2112) 


0- 

1- 


'So' 

's. 


IS (2076) 


1 
2 


-1 


2072 


-4 


Dsi (2535) 
Ds2 (2573) 


1+ 
2+ 


'Pi ' 

'P2 


IP (2559) 


3 
2 


-2 


2563 


4 


bu, bd quarks 
















B (5279) 
B* (5325) 


0" 
1- 


'So 
'S, _ 


IS (5314) 


1 
2 


-1 


5315 


1 


bs quarks 
















Bs (5374) 
B* (5421) 


0" 

1- 


'So' 
'S, 


IS (5409) 


1 
2 


-1 


5408 


-1 
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FIGURES 



Figure 1: Dependence on the basis state scale (3 of the three lowest positive energy 
states and the three highest negative energy states for the no-pair equation with 
flux tube. Here we have chosen niq = 0.3 GeV, a = 0.2 GeV^, k, = 0.5, k = —1 and 
j = The full lines correspond to iV = 25, and dashed lines correspond to iV = 15 
basis states used. 

Figure 2: Regge trajectories for the no-pair equation with flux tube. Again, we have 
chosen rriu^a = 0.3 GeV, a = 0.2GeV^, and n = 0.5. The full lines correspond to 
/c = — (j -|- |), and dashed lines to = j + |. To ensure that all calculated energies 
are correct, we used = 50 basis states, and kept first 15 states. 

Figure 3: The Isgur-Wise function for B D^*^li'i decays calculated from the no- 
pair equation with tube confinement. Values for the light quark mass rriu^d, tension 
a and short range potential constant k are taken from the fit (|65D. For the sake of 
clarity, error bars are shown only for the CLEO II data. 
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Figure 1 
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Figure 2 



20 




Figure 3 
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